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ABSTRACT 

We present a spectroscopic catalog of 70,841 visually inspected M dwarfs from the seventh data 
release of the Sloan Digital Sky Survey (SDSS). For each spectrum, we provide measurements of the 
spectral type, a number of molecular bandheads, and the Ha, H/3, H7, HJ and Ca II K emission lines. 
In addition, we calculate the metallicity-sensitive parameter Q and identify a relationship between Q 
and the g — r and r — z colors of M dwarfs. We assess the precision of our spectral types (which 
were assigned by individual examination), review the bulk attributes of the sample, and examine 
the magnetic activity properties of M dwarfs, in particular those traced by the higher order Balmer 
transitions. Our catalog is cross-matched to Two Micron All Sky Survey (2MASS) infrared data, and 
contains photometric distances for each star. Lastly, we identify eight new late-type M dwarfs that are 
possibly within 25 pc of the Sun. Future studies will use these data to thoroughly examine magnetic 
activity and kinematics in late-type M dwarfs and examine the chemical and dynamical history of the 
local Milky Way. 

Subject headings: stars: low-mass — brown dwarfs — stars: activity — stars: late-type — stars: 
abundances — Galaxy: kinematics and dynamics 



1. INTRODUCTION 

Over the past decade, wide-field, deep astronomical 
surveys have provided an unprecedented statistical plat- 
form for studying the Universe (e.g. Sloan Digital Sky 
Survey, Two Micron All Sky Survey, 2dF, UKIRT In- 
frared Deep Sky Survey). For the lowest- mass and most 
populous stars in the Milky Way (M dwarfs), the Sloan 
Digital Sky Survey (SDSS) has yielded photometric sam- 
ples that exceed 30 million stars (e.g. Bochanski et al. 
2010) and spectroscopic samples of over 40,000 M and L 
dwarfs (West et al. 2008). Previous SDSS studies have 
elucidated the mean properties of low-mass stars (Haw- 
ley et al. 2002; West et al. 2004, 2005; Bochanski et al. 
2007b; Schmidt et al. 2010b), their magnetic activity and 
flaring properties as a function of mass (West et al. 2004, 
2008; Kowalski et al. 2009; Kruse et al. 2010; Hihon et al. 
2010), as well as the structure and kinematics of the local 
Milky Way thin and thick disks (Bochanski et al. 2007a; 
Juric et al. 2008; Fuchs et al. 2009; Bochanski et al. 2010) 
and the low-mass initial mass and luminosity functions 
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(Covey et al. 2008b; Bochanski et al. 2010). 

M dwarfs have main sequence lifetimes that are con- 
siderably longer than the age of the Universe and can be 
used to trace the evolution of both stellar properties and 
the Milky Way disks. West et al. (2006, 2008) showed 
that magnetic activity (as traced by Ha) in M dwarfs 
decreases with age and that M dwarfs appear to have 
finite activity lifetimes from '^1-2 Gyr for early- type M 
dwarfs (M0-M3) to -7-8 Gyr for later type stars (M5- 
M7). These results are important for the habitability of 
extrasolar planets orbiting M dwarfs (e.g. Charbonneau 
et al. 2009), as active stars may disrupt planetary at- 
mospheres (Segura et al. 2010). The Ha emission does 
not have sufficient energy to significantly affect planetary 
systems but has been found to be correlate with X-ray 
emission (Reid et al. 1995b; Covey et al. 2008a), which 
can interact with extrasolar planet atmospheres. 

While Ha is the most commonly studied emission line 
in M dwarfs, the higher-energy hydrogen Balmer and Ca 
II transitions are also present in the optical spectra of 
active stars. The higher energy emission lines appear to 
trace different temperature regions in the chromosphere 
and can be used to characterize the upper atmospheres 
of active M dwarfs (Walkowicz & Hawley 2009, and ref- 
erences therein). However, lacking large samples of low- 
mass stars with spectroscopic coverage across the entire 
optical bandpass, it is still not clear how the various 
activity indicators trace each other in active M dwarfs 
(Rauscher & Marcy 2006; Walkowicz & Hawley 2009). 

West et al. (2008, hereafter W08) also showed that the 
ratio of CaH/TiO molecular indices (a quantity that is 
likely related to metallicity; Gizis 1997) decreased with 
height in the Galactic disk (a proxy for age). This tan- 
talizing result suggests that M dwarfs can be used to re- 
construct the chemical evolution of the local Milky Way 
disk. However, the inability to assign accurate metallici- 
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ties to individual M dwarfs limits a detailed exploration 
of Galactic chemical evolution. Several recent studies 
have made progress in calibrating the M dwarf metallic- 
ity scale. Rojas-Ayala et al. (2010) used infrared spec- 
troscopy of M dwarfs in wide binaries with higher-mass 
stars to calibrate infrared spectral features with metal- 
licity. In the optical, Lcpine et al. (2007) developed a 
metallicity-dependent quantity (, which is a relation be- 
tween the CaH and TiO molecular indices and defined 



as: 



where 



c = 



1 - Ti05 
- [Ti05]z, 



(1) 



[Ti05]2 



-0.164(CaH2-FCaH3)3 

+0.670(CaH2 + CaHS)^ 
-0.118(CaH2 + CaH3) - 0.050. 



(2) 



Using wide binaries consisting of an M dwarf and a higher 
mass star, Woolf et al. (2009) showed that the ( index 
is a good discriminant between high and low-metallicity 
([Fe/H] =0 and [Fe/H] = -1 respectively) for early- type 
M dwarfs (-M0-M3). While ( is a useful tool for 'ex- 
amining the relative metallicity among stars, finer deter- 
minations of an absolute metallicity from C will require 
additional calibration. 

Before the SDSS (York et al. 2000), the largest spectro- 
scopic samples of M dwarfs contained only a few thou- 
sand individual stars (Reid et al. 1995a) and were pri- 
marily focused on the red portion of the optical spec- 
trum (~6000-8000 A). The spectroscopic catalogs of M 
dwarfs from the first few SDSS data releases quickly sur- 
passed previous records for sample size (Hawley et al. 
2002; West et al. 2004) and had spectral coverage that 
spanned the entire optical band (~3900-9300 A) at low 
resolution (R~1800). W08 presented a spectroscopic cat- 
alog of more than 44,000 M dwarfs from the SDSS Data 
Release 5 (DR5), more than doubling the M dwarf spec- 
troscopic tally. Due to the large number of spectra in 
the DR5 sample, spectral types were assigned using the 
Hammer^ ^automatic spectral typing facility (Covey et al. 
2007). While the Hammer spectral types are generally 
good to within ±1 sub type, there have been recent in- 
dications that there may be systematic offsets for some 
of the M dwarf spectral types (see Section 2). There- 
fore, to ensure the highest quality sample, it is impor- 
tant to visually inspect each spectrum. There have been 
two subsequent SDSS data releases since the W08 sample 
that include over 50,000 additional M dwarf candidates 
and contain many sightlines at low Galactic latitudes (as 
part of the SEGUE survey; Yanny et al. 2009). The new 
lines of sight extend both the radial and vertical extent 
of the M dwarf sample, and provide a larger statistical 
platform from which to probe the structure, kinematics 
and evolution of the Milky Way using its smallest stellar 
constituents. 

Here we present the latest spectroscopic catalog of M 
dwarfs from Data Release 7 (DR7; Abazajian et al. 2009) 
of the SDSS. While the main thrust of this paper is the 

The Hammer software is available for public download at: 
http: / / astro.washington.edu/users/slh/hammer 



presentation and characterization of a new spectroscopic 
sample, we have included a general magnetic activity 
analysis as well as a discussion about the future use of 
this sample to probe Galactic chemical evolution. We 
describe our sample selection techniques in Section 2, 
specifically discussing our manual inspection of all M 
dwarf candidates. Section 3 gives the mean properties 
of the sample, including a new activity analysis of the 
hydrogen Balmer and Ca II K emission features. Section 
3 also contains a list of candidate M dwarfs within 25 
pc and explores how the metallicity sensitive parameter 
C varies as a function of stellar parameters. We discuss 
our results and future use of our new sample in Section 
4. 

2. DATA AND SAMPLE SELECTION 

The SDSS provides a large, uniform photometric and 
spectroscopic dataset from which to extract high-quality 
samples of low-mass dwarfs (Gunn et al. 1998; Fukugita 
et al. 1996; Hogg et al. 2001; Gunn et al. 2006; Ivezic 
et al. 2004; Pier et al. 2003; Smith et al. 2002; Tucker 
et al. 2006). We used the DR7 CasJobs^^ ^q^j ^^i^^^^ 
116,161 M dwarf candidate objects with SDSS spectra 
that 1) fell within the typical color range for M and L 
dwarfs (r-i > 0.42 and i-z > 0.24; W08, Kowalski et al. 
2009); 2) had SDSS spectral classifications of "STAR" 
or "STARXATE" ; and 3) had radial velocities smaller 
than 1500 km s'^ Kowalski et al. (2009) found that 
the color cuts used in previous studies did not extend 
blue enough (r — « in particular) to include all of the MO 
dwarfs. Our new (bluer) color selection corrects for this 
previous oversight. 

The spectra for all 116,161 candidates were visually in- 
spected and assigned spectral types using the Hammer 
spectral typing facility (Covey et al. 2007). Low signal- 
to-noise ratio (SNR) spectra (SNR <3 at ~8300 A) and 
extragalactic interlopers were removed during visual in- 
spection, resulting in 109,639 objects. We also removed 
stellar spectra that were not identified as M dwarfs dur- 
ing visual inspection as well as those spectra that were 
duplicated in our sample (including 6771 spectra for 2661 
M dwarfs that were taken on different days) and arrived 
at the 70,841 spectra of M dwarfs in our DR7 catalog. 
While many of the objects that we removed were K or 
L dwarfs, others were white dwarf M dwarf (WD-dM) 
pairs and low-metallicity subdwarfs. All of the L dwarfs 
that were removed from our DR7 sample were cataloged 
and analyzed by Schmidt et al. (2010b). Many of the 
other objects that were not included in our DR7 sample 
will be presented in future studies. Not all of the WD- 
dM pairs were removed during the visual inspection. We 
used the color cuts from Smolcic et al. (2004, u — g < 
2, g — r > 0.3, r — i > 0.7, cTu,g,r,i < 0.1) to identify 
additional WD-dM pairs (that avoided manual detec- 
tion). All of the 497 objects that matched these color 
criteria are flagged as "WDM" in our DR7 catalog. We 
also defined a clean photometric sample using the SDSS 
photometric processing flags (SATURATED, PEAK- 
CENTER, NOTCHECKED, PSF.FLUX.INTERP, IN- 
TERP .CENTER, BAD.COUNTS.ERROR were afl set 
to zero in the r, i, and z bands)^'', which resulted in 

http://casjobs.sdss.org/ 

http://www.sdss.org/dr7/products/catalogs/flags.html 



SDSS DR7 M Dwarfs 



3 



65,277 M dwarfs. While our catalog contains all of the 
70,84 visually inspected M dwarf spectra, the "GOOD- 
PHOT" and "WDM" flags can be used to obtain samples 
that include good photometry and remove possible WD- 
dM pairs respectively. 

Radial velocities (RVs) were measured by cross- 
correlating each spectrum with the appropriate Bochan- 
ski et al. (2007b) M dwarf template. This method has 
been shown to produce uncertainties ranging from 7-10 
km s-i (Bochanski et al. 2007b). All of the DR7 objects 
were cross-matched to the USNO-B/SDSS proper mo- 
tion catalog (Afunn et al. 2004, 2008), identifying 39,151 
M dwarfs with good proper motions^'*. Distances to 
each star were calculated using the Mr vs. r — z color- 
magnitude relation given in Bochanski et al. (2010). Our 
calculated distances have uncertainties of ~20%, arising 
mostly from the intrinsic spread of the main sequence. 
The proper motions and distances were combined with 
the RVs to produce 3-dimerisionaI space motions for the 
DR7 M dwarfs. Although we include the standard U, 
V, W space motions in our catalog, we caution that the 
U, V, W velocities are in a Cartesian coordinate sys- 
tem that may not be appropriate for stars at apprecia- 
ble distances from the Sun. We therefore also include 
the Galactic radial {R), tangential (6) and vertical {Z) 
cylindrical components of the position and velocity for 
each star for which we have 3D space motions. 

We also matched our catalog to the 2MASS point 
source catalog (Cutri et al. 2003), matching only to 
unique 2MASS counterparts within 5" of the SDSS po- 
sition that do not fall within the boundaries of an ex- 
tended source (gal_CONTAM =0). To ensure that we 
used only high quality 2MASS photometric data, we 
applied additional cuts to each of the J, and Ks 
bands. If the source was not detected (rd_flg =0), 
nominally detected (rd_flg =6), was detected but unre- 
solved (rd_flg =9), or had contaminated/confused pho- 
tometry (cc_FLAG 7^0) in a particular band, the 2MASS 
data were not included. This resulted in 57,956 2MASS 
counterparts with J, H and Kg magnitudes and their 
uncertainties that were included in the catalog. 

We used the 2MASS photometry to investigate any 
possible M giant contamination. Bessell & Brett (1988) 
found that M giants and M dwarfs separate in J — H vs. 
H — K color space due to differences in H2O absorption 
in their atmospheres. By comparing the DR7 J — H 
vs. H — K color-color diagram with that of Bessell & 
Brett (1988), we find that no more than 0.5% of our 
sample could be giants. In addition. Covey et al. (2008b) 
conducted a complete magnitude-limited spectroscopic 
survey of a in° field in the SDSS footprint and concluded 
that the giant contamination rate was less than 2% for 
stars redder than a spectral type of K5. 

2.1. Spectral Typing 

Previous large spectroscopic samples of M dwarfs have 
relied on automated spectral typing due to the quantity 
of time required to manually inspect tens of thousands 
of candidate spectra. Systematic discrepancies were re- 
cently identified in the spectral types from the W08 

Good proper motions arc defined as those with MATCH =1, 
DIST22 > 7, SIGRA < 1000, SIGDEC < 1000 and (NFIT = 6 or 
(NFIT = 5 and (O < 2 or J < 2))) (Munn et al. 2004). 
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Fig. 1. — The distribution of differences between the Hammer 
automatic spectral types and the mean of the visually inspected 
spectral types for M dwarfs with visual types earlier than M5 (left) 
and M5 and later (right). While there is good agreement for the 
early-type M dwarfs, there is a systematic offset in the automatic 
types assigned by the Hammer; the Hammer assigns an earlier type 
~38% of the time. 



SDSS sample that were determined automatically by the 
"Hammer" spectral typing algorithm (Covey et al. 2007). 
The bias was detected as a systematic offset for late-type 
stars, whose automatic classification was often 1 subtype 
earlier than determined via visual inspection (see Fig- 
ure 1). We thus decided that visual inspection would 
produce the most reliable and precise spectroscopic sam- 
ple. We visually inspected all 116,161 M dwarf candi- 
dates and manually assigned spectral types. The sample 
was divided among 17 individuals^^ who used the man- 
ual "eyechcck" mode of the Hammer (v. 1_2_5) to assign 
spectral types and remove non-M dwarf interlopers. 

Figure 1 shows the difference between the spectral 
types automatically determined by the Hammer and the 
mean visual inspection. For early-type M dwarfs (left 
panel), the Hammer and the spectral typcrs agree most 
of the time. However, for -^38% of the late- type M dwarfs 
(right panel), the Hammer assigns spectral types 1 sub- 
type earlier than the average human spectral typer. This 
confirms that while the Hammer generates automatic 
spectral types within the quoted ±1 subtype accuracy 
(Covey et al. 2007), there is a systematic offset of 1 sub- 
type for the late-type stars, and justifies our effort to 
manually inspect more than 116,000 candidates. 

Each spectral typer also examined a control sample of 
1000 spectra, of which 638 were M dwarfs. We used this 
control sample to assess the quality and reliability of the 
visually inspected sample by quantifying the variations 
among the 17 individual typers, and by comparing the 
median visual type for each star to its automatic Ham- 
mer type. The results of the control sample show that 
there is excellent agreement between all of the spectral 
typers with a large fraction of the stars being assigned. 
Almost all of the visual classifications in the control sam- 
ple agreed to within ±1 subtype of the median value; for 
most of the stars, the dispersion in visual classifications 
was < 0.4 subtypes. To ensure that there were no major 
systematics in the spectral types assigned by any indi- 
vidual, we examined the mean difference between each 
typer's visual classifications of the control sample spec- 
tra and the median classifications of the entire group. 
We divided the sample into the early types (< M5) and 

The order of the co-authors was based on the number of spec- 
tra examined. 
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late types (> M5) to investigate any spectral type de- 
pendence in typing quality: all of the spectral typers 
clustered near the median values with no individual dis- 
playing a systematic offset larger than 0.2 subtypes. 

2.2. Measured Quantities 

As part of our analysis wc; mcasiircd a number of spec- 
tral lines and molecular features in each M dwarf spec- 
trum. All of the spectral measurements were made using 
the RV corrected spectra. The TiOl, Ti02, Ti03, Ti04, 
Ti05, Ti08, CaOH, CaHl, CaH2, and CaH3 molecu- 
lar bandhead indices and their formal uncertainties were 
measured using the Hammer with the molecular band- 
heads as defined by Rcid et al. (1995a) and Gizis (1997). 
We also measured the chromospheric hydrogen Balmer 
and Ca II lines that are associated with magnetic activ- 
ity. We expanded the Ha analysis of West et al. (2004, 
2008) to include H/3, H7, R6 and Ca II K (He and Ca II 
H are blended in SDSS data and were not included in our 
sample) . All of the line measurements were made by in- 
tegrating over the specific line region (8 A wide centered 
on the line) and subtracting off the mean flux calculated 
from two adjacent continuum regions (Continuum A and 
Continuum B in Table 1). Equivalent widths (EW) were 
computed for each line by dividing the integrated line 
flux by the mean continuum value (as in previous stud- 
ies, we define the EW to be positive for emission lines; 
Hawley et al. 1996; West et al. 2004). The low-resolution 
of SDSS spectra does not affect our ability to measure 
accurate line values, with the exception of Call K. Due to 
the intrinsic Call K absorption present in M dwarf photo- 
spheres, weak activity appears as an emission peak in the 
absorption line. This weak emission cannot be resolved 
in SDSS spectra. However strong Call K emission over- 
whelms the absorption and can be easily recovered from 
low-resolution spectra (Walkowicz & Hawley 2009). The 
lines and continuum regions (in vacuum wavelengths) can 
be found in Table 1. Formal uncertainties on the EWs 
for each line were also computed. 

Similar to the West et al. (2004, 2008) samples, we 
defined magnetically active stars as those that had de- 
tectable emission lines in their spectra. Previously, this 
had been done solely for Ha. To be classified as active in 
a specific emission line, the following four criteria had to 
be satisfied: 1) the EW of the line was larger than some 
minimum value (see below); 2) the EW value must be 3 
times the uncertainty; 3) the signal to noise ratio (SNR) 
in the continuum must be larger than 3; and 4) the height 
of the spectral line must be larger than 3 times the noise 
in the continuum. Stars were classified as inactive if they 
met criterion #3 and had no detectable emission. Crite- 
rion =f/=3 preferentially selects brighter stars at any given 
color or spectral type. This selection bias should not have 
a large effect our results since we examine the properties 
of our sample in a Galactic context and both active and 
inactive stars with low SNR should be removed. How- 
ever, if active stars are intrinsically more luminous than 
their inactive counterparts (at the same color; see Sec- 
tion 3.2 below), criterion ^3 may serve to include slightly 
more active stars in our resulting analyses. 

Because the strength of the stellar continuum changes 
dramatically from Call K to Ha in M dwarfs, the EW 
thresholds used in criterion ^1 are different for each 
line. We experimentally derived the EW thresholds by 



automatically measuring activity in a subsample of 300 
M dwarfs (that spanned a range of SNR and spectral 
types) using 10 different EW activity thresholds. Addi- 
tionally, each spectrum was visually inspected and the 
activity state for all five emission lines was determined. 
By comparing the results of our manual inspection with 
the 10 different automatic routines, we determined the 
EW thresholds that most accurately reproduced our vi- 
sual activity classifications. The EW activity threshold 
for each of the five emission lines can be found in Table 
1. We found that the previous Ha activity limit of 1 A 
(West et al. 2004) was slightly too high and excluded a 
number of good detections. We have therefore reduced 
the Ha EW activity criterion to 0.75 A. 

Recently, Kruse et al. (2010, hereafter KIO) investi- 
gated the time variability of SDSS M dwarfs and defined 
their Ha activity thresholds based on the median EW 
values that were detected at a 3a confidence level as a 
function of spectral type. Although the KIO threshold 
EWs provide a robust way to select a clean sample, it 
unduly biases the sample toward the most active stars 
at later spectral types. For example, our DR7 catalog 
contains a large number of active M6 dwarfs (that have 
Ha emission lines with peak values that are more than 

5 times larger than the noise) that have EWs below the 
3A threshold used in the KIO study. In addition, the 
activity fractions presented in KIO (see their Figure 4) 
include these active stars (with EW < 3A) in the "in- 
active" component of the activity fraction denominator. 
This explains why the activity fractions reported by KIO 
are smaller than previous determinations (West et al. 
2004, 2008) and indicates a serious selection bias in the 
KIO results. 

For all of the active stars in the sample we computed 

the ratio of luminosity in the emission line as compared to 
the bolomctric luminosity (Lunc/iboi)- Wc followed the 
methods of Hall (1996), Walkowicz et al. (2004) and West 

6 Hawley (2008) who derived x factors for the Hydrogen 
Balmer and Ca II chromospheric lines as a function of M 
dwarf spectral type. The x factor uses empirical results 
for the bolomctric luminosities of M dwarfs as a function 
of color and/or spectral type (Leggett et al. 1996, 2001) 
and relates the EW of a line to the fraction of the bolo- 
metric luminosity emitted by the line. The iune/iboi 
values were computed by multiplying the EW of each 
active star by the appropriate x value. Formal uncer- 
tainties were computed for each Lune/ -^boi value and are 
included in the final database. 

We also computed the metal sensitive parameter de- 
fined by Lepine et al. (2007), which uses a combination 
of the Ti05, CaH2, and CaH3 molecular band indices to 
separate the sample into different metallicity classes (see 
Equations 1 & 2). This is similar to the Gizis (1997) clas- 
sification system but was re-calibrated using wide com- 
mon proper motion pairs that were assumed to be at the 
same metallicity. Stars with solar metallicity ([Fe/H]=0) 
have C values ~1 and stars with [Fe/H]=-1 have ^OA 
(Woolf et al. 2009). Although there is considerable scat- 
ter in the [Fe/H] versus ( relation at high-metallicities, 
this parameter is very useful for finding and classify- 
ing low-metallicity stars that are likely members of the 
Galactic halo. 

Our catalog containing all of the measured quantities is 
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TABLE 1 
Emission Lines 



Line 


Central 


Continuum A 


Continuum B 


Activity Limit 




Wavelength (A) 


(A) 


(A) 


EW (A) 


Ha 


6564.66 


6555.0-6560.0 


6570.0-6575.0 


0.75 


H/3 


4862.69 


4840.0-4850.0 


4875.0-4885.0 


1.00 


H7 


4341.69 


4310.0-4330.0 


4350.0-4370.0 


1.00 


m 


4102.90 


4075.0-4095.0 


4110.0-4130.0 


1.50 


Call K 


3934.78 


3952.8-3956.0 


3974.0-3976.0 


1.50 



Note. — All of the wavelengths are given in vacuum units 



publicly available on the Vizier site or can be obtained 
by contacting the corresponding author. The individual 
spectra are available from the SDSS DR7 website^''. As 
with previous SDSS spectroscopic catalogs of low-mass 
stars, we remind the community that these data do not 
represent a complete sample and that the complicated 
SDSS spectral targeting (see Section 3.1) introduces a va- 
riety of selection effects. However, our new sample covers 
a large range of values for many of the physical attributes 
of the M dwarfs, including parameters that arc sensitive 
to activity, metallicity, and Galactic motion, making ac- 
curate activity, kinematic, and chemical analyses possi- 
ble. In addition, because some of the derived quantities 
are computed by automatic routines, values for a small 
percentage (^4%; West ct al. 2004) of individual stars 
may be incorrect; this should not affect large statistical 
results. Users are nevertheless cautioned to understand 
the origin of specific data products before using them 
indiscriminately. 

3. RESULTS 

Figure 2 shows the distribution of spectral types for the 

DR7 low-mass star sample. We included the entire visu- 
ally confirmed sample (solid), the stars for which we have 
good photometry (as indicated by our GOODPHOT flag; 
dotted), and the stars for which we have 3D kinematic 
information (as set by our GOODPM flag; dashed). A 
majority of the late-type M dwarfs are removed from the 
kinematic sample because they are too faint to be de- 
tected in the USNO-B catalog. Future studies will mea- 
sure proper motions for the late-type M dwarfs in the 
DR7 catalog and create a more representative kinematic 
sample. 

3.1. Median Colors 

All of the magnitudes included in our DR7 catalog and 
used to compute the median colors were extinction cor- 
rected using the Schlegel et al. (1998) linc-of-sight dust 
maps (and the dust model of Cardelli et al. 1989, and Ry 
= 3.1) . We provide extinction corrections for each band 
in the catalog so that observed magnitudes can be eas- 
ily reproduced. Because of the close proximity of many 
of the M dwarfs in the sample, the entire Schlegel ct al. 
(1998) extinction is likely an overestimate for many lines- 
of-sight. The uncertainty in the extinction correction 
introduces a small uncertainty in the derived distances. 
However, given the small amount of extinction toward 

http;//vizier.u-strasbg.fr/cgi-bin/VizieR 
^"^ http://www.sdss.org/dr7 
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Fig. 2. — The distribution of spectral types in the spectroscopic 
SDSS DR7 M dwarf catalog (solid). The dotted line represents the 
stars for which good quality photometry exists (GOODPHOT = 
1) and the dashed line indicates the stars for which good proper 
motions exist (GOODPM =1). Due to the shallow depths and blue 
sensitivity of the USNO-B catalog, there are only a small number 
of measured proper motions of the latest-type M dwarfs. 



objects in the SDSS footprint, the added uncertainties 
to the distances are typically less than 5%. We are cur- 
rently determining the amount of extinction needed to 
reproduce an unreddened spectrum for each object in 
our DR7 sample, and will report on that analysis in a 
future paper (Jones et al. 2011). 

We computed median colors for each subtype of the 
DR7 sample, applying stringent cuts to include only 
those stars with the highest quality photometry. To be 
included in the median color calculation, SDSS stars had 
to have GOODPHOT = 1, r-band extinction <0.05 mag- 
nitudes, and r, i and 2;-band magnitude uncertainties 
<0.05 magnitudes. Median colors incorporating 2MASS 
magnitudes were calculated from the subset of stars with 
J, H, and Ks uncertainties <0.05. The median colors 
and standard deviations (in parentheses) are shown in 
Table 2. Due to the systematic offset in the Hammer 
automatic spectral types for late-type M dwarfs, our av- 
erage colors are bluer than previously reported (W08, 
Kowalski ct al. 2009). We include the g — r colors but 
warn that they are likely metallicity dependent (West 
et al. 2004; Bochanski et al. 2007b, see also Section 3.3 
below). 

Figure 3 shows the r versus r — z Hess diagram for 
the M dwarfs in the SDSS DR7 sample. The gaps in 
color space reflect the non-uniform sampling due to the 
spectral targeting algorithm. To quantify the effect of 
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TABLE 2 
Median Colors 



Sp. Type AfsDSs'' iV2MASs'' 9-r'' r-i i-z z-J J -H H - Ks 



MO 915 259 1.31 (0.16) 0.56 (0.08) 0.33 (0.06) 1.20 (0.25) 0.64 (0.07) 0.16 (0.07) 

Ml 699 230 1.39 (0.16) 0.73 (0.10) 0.41 (0.08) 1.32 (0.41) 0.62 (0.08) 0.19 (0.08) 

M2 1078 689 1.40 (0.13) 0.96 (0.10) 0.53 (0.08) 1.26 (0.18) 0.60 (0.07) 0.22 (0.08) 

M3 1220 831 1.41 (0.13) 1.13 (0.11) 0.61 (0.08) 1.30 (0.24) 0.59 (0.06) 0.23 (0.06) 

M4 814 522 1.46 (0.15) 1.33 (0.13) 0.71 (0.09) 1.37 (0.15) 0.59 (0.06) 0.25 (0.06) 

M5 328 269 1.52 (0.13) 1.62 (0.21) 0.90 (0.13) 1.47 (0.12) 0.59 (0.06) 0.28 (0.07) 

M6 359 401 1.56 (0.15) 1.92 (0.14) 1.05 (0.07) 1.61 (0.08) 0.60 (0.06) 0.31 (0.08) 

M7 160 352 1.58 (0.16) 2.09 (0.19) 1.14 (0.10) 1.71 (0.10) 0.60 (0.08) 0.34 (0.07) 

M8 6 118 1.64 (0.12) 2.56 (0.24) 1.41 (0.13) 1.93 (0.13) 0.63 (0.07) 0.39 (0.08) 

M9 4 77 1.87 (0.74) 2.70 (0.20) 1.71 (0.18) 2.15 (0.14) 0.69 (0.09) 0.42 (0.07) 



Note. — The median and rms scatter (in parcnilicscs) of the extinction corrected SDSS-2MASS colors for each M dwarf 



spectral type, u — g colors were not included due to the large uncertainties in the u-band photometry for most of the sample 
M dwarfs. The g — r colors are included but are likely metallicity dependent (West et al. 2004; Bochanski et al. 2007b). 

^The SDSS numbers reflect the stars that have good SDSS photometry (the GOODPHOT flag set to 1) and have r, i and z 
photometric uncertainties < 0.05 magnitudes. 

^The 2MASS numbers indicate the number of stars that have z, J, H and its-band uncertainties < 0.05 magnitudes. 

*^The g — r colors were calculated using stars with good SDSS photometry (the GOODPHOT flag set to 1) and g and r-band 
photometric uncertainties < 0.05 magnitudes. 



SDSS spectroscopic targeting algorithms on the median 
colors, we conducted the following test. We began with 
the straw man assumption that the median r — i, i — z 
and r — z color-spectral type relations reported in Table 
2 are representative of the underlying stellar population. 
We then used these quantities, along with the observed 
color spread in each spectral type bin, to construct a 
synthetic stellar population with a uniform distribution 
in both magnitude and spectral type. We then applied 
photometric cuts to the smooth underlying distribution, 
replicating the SDSS spectroscopic selection and produc- 
ing a synthetic "observed" sample that is highly struc- 
tured in color-magnitude space. The median colors and 
spreads were then calculated as a function of spectral 
type from this synthetic sample and compared to the 
measured colors in Table 2. In general, the colors of the 
synthetic stars agree to within 0.02 mags of the observed 
sample (an agreement that would not be expected if the 
underlying population were different from the observed 
SDSS stars ) and we conclude that the colors presented 
in Table 2 are not severely affected by the spectroscopic 
selection. 

3.2. Activity 

Figure 4 shows the Ha activity fraction as a function 
of M dwarf spectral type for stars with data quality suf- 
ficient to measure Ha emission (and GOODPHOT =1 
and WDM =0). Our results are in good agreement with 
previous studies (West et al. 2004, W08), but with much 
lower uncertainties due to the ^59,000 M dwarfs used 
to generate Figure 4 (7952 are Ha active and 51,034 arc 
Ha inactive). There are a few small differences between 
Figure 4 and the Ha activity fractions reported by W08 
that are likely due to the changes in some of the spectral 
types (see above) and the different Galactic sightlines in- 
cluded in the DR7 sample. Because of the lower SNR in 
the blue portion of M dwarfs, fcnvcir stars were available 
to measure activity for the higher order Balmer and CaH 
K emission lines. Table 3 gives the number of active and 
inactive stars for each emission line indicator that had 
the necessary quality for our activity analysis. 

As discussed in W08, the activity fraction is highly de- 




r-z 



Fig. 3. — The r versus r — z Hess diagram for the M dwarfs in the 
SDSS DR7 sample. The gaps in color space reflect the iion-uiiiforin 
sampling due to the spectral targeting algorithm. The irregular 
sampling has a minimal effect on the median colors derived in the 
DR7 sample. 



pendent on the location of the samples in the Galaxy. 
Because M dwarfs have finite activity lifetimes and are 
dynamically heated in the Galactic disk as they age, their 
activity state is correlated with position in the Galaxy. 
One of the reasons that the M dwarfs in our sample have 
much smaller activity fractions than those studied nearby 
(e.g. Hawley et al. 1996; Gizis et al. 2000), is that the 
SDSS volume is much larger than those used in previous 
catalogs and is concentrated on the north galactic cap, 
yielding a much older stellar population. This is par- 
ticularly clear for the early-type M dwarfs, which have 
median distances greater than 500 pc, and therefore ages 
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Fig. 4. The Ha magnetic activity fraction as a function of spec- 
tral type. Error bars were calculated from the binomial distribu- 
tion and the number above each data point indicates the numbers 
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TABLE 3 

ACTIMT^' iNniCATORS 



Line 


-^act i vc 


-^inactive 


Active Fraction 


Ha 


7952 


51034 


0.13 


H^ 


2236 


35144 


0.06 


H7 


1175 


18816 


0.06 


H5 


528 


14664 


0.03 


CaH K 


620 


9053 


0.06 



that arc considerably older than their short activity life- 
times (~ 1-2 Gyr). We reiterate the warning from W08 
that activity fractions in M dwarfs must be discussed in 
the proper Galactic context. 

The various emission lines measured in our spectra 
arc formed at slightly different locations in the chro- 
mosphere, suggesting that the strength of one emission 
line may not necessarily predict the strength of another. 
We therefore examined how the various activity induced 
emission lines trace each other as a function of spectral 
type and absolute distance from the Galactic plane (a 
proxy for age). Figure 5 shows the Ha (diamonds) and 
H/3 (asterisks) activity fractions for M2-M7 dwarfs as a 
function of Galactic height. Only the Ha activity frac- 
tions for stars that could have detected H/3 emission are 
included. This limits the sample to the brighter stars at 
each spectral type that meet the SNR criterion and does 
not include all of the M dwarfs that are sensitive to Ha 
emission. Spectral types earlier than M2 and later than 
M7 were not included due to an insufficient sample size 
of active stars. As seen in previous studies, the mag- 
netic activity fractions decrease as a function of Galac- 
tic height, confirming an age- activity relationship; stars 
closer to the Plane are statistically younger and more 
likely to be active. Figure 5 also demonstrates that Ha 
and H/3 trace each other extremely well for the entire 
DR7 sample, implying that the H/3 activity lifetimes are 
essentially the same as those for Ha. 

The other activity tracers also appear to correlate with 
Ha emission. Figure 6 shows the activity fractions for 
all five emission lines in M4 dwarfs as a function of ver- 
tical distance from the Galactic plane. Each panel plots 
the activity fractions for stars where the data are good 



enough to measure (i.e. are sensitive to) H/3 (top left), 
H7 (top right), US (bottom left) and CaH K (bottom 
right) activity as well as the emission lines redder than 
that specific tracer. The activity fractions for all of the 
lines are in excellent agreement with each other and indi- 
cate that all five emission lines are produced for the same 
amount of time during the lifetime of an M dwarf. This 
new result is complementary to previous studies that 
found strong correlations among different emission line 
strengths in M dwarfs (Rauschcr & Marcy 2006; Walkow- 
icz & Hawley 2009). Although the relative emission line 
strengths for individual stars may fluctuate over time (see 
Cincunegui et al. 2007), the mean values of the emission 
lines are well correlated in the large DR7 sample. The 
unprecedented size and Galactic distribution of our spec- 
troscopic sample confirms that not only the strength, but 
the duration of activity appears to be similar for all of 
the Balmer and CaH activity indicators. 

Figure 6 also demonstrates a possible (and important) 
selection effect present in magnitude-limited or activity- 
selected samples. A close examination of Figure 6 (and 
Figure 5) reveals that the activity fractions arc larger for 
the stars that are sensitive to the bluer emission lines. 
This trend was not expected since the stars used to com- 
pute the activity fraction were drawn from the same vol- 
umes in each of the four panels. The only difference was 
that the stars selected for the bluer emission line analy- 
ses were required to have higher SNR. With this in mind, 
we can explain the different activity fractions as follows. 

The actual spread in absolute magnitude (or luminos- 
ity) at a given spectral type can be quite significant.^^ 
Some of this spread may be due to differences in the phys- 
ical properties of the stars. Hawley et al. (1996) showed 
that magnetically active stars are brighter in My than 
their inactive counterparts. Recently, Bochanski et al. 
(2011, hereafter Paper H) found that both active and 
higher metallicity M dwarfs (many stars are both) appear 
to be brighter in Mr at a given color or spectral type. Be- 
cause the bluer emission lines that wc; are examining in 
Figure 6 require higher SNR spectra to accurately mea- 
sure activity, the M4 dwarfs at a given distance that 
meet the SNR criterion are preferentially more luminous 
than other stars in the same distance bin. If activity is 
correlated with luminosity, our SNR selection will bias 
our sample toward containing a higher fraction of active 
stars. 

In addition to activity fractions, we also investigated 
the median luminosity in each emission line as a func- 
tion of spectral type. Previous studies have investigated 
how the Ha luminosity changes as a function of spectral 
type (Hawley et al. 1996; Gizis et al. 2000; Burgasser 
et al. 2002; West et al. 2004). These studies found that 
the fraction of luminosity emitted in Ha is constant over 
the range M0-M5 and decreases by almost an order of 
magnitude for late- type M dwarfs. However, most of the 
previous studies concentrated on Ha and did not include 
similar analyses for the other Balmer or Call lines, and 
those that did were based on a small sample of stars 
(Hawley & Pettersen 1991; West & Hawley 2008). 

Recent studies have confirmed that the spread in absolute 
magnitude as a function of color is much smaller (Bochanski ct al. 
2010, 2011; Fahcrty ct al. 2011). This is the primary reason that 
our distances are estimated from colors and not spectral types. 
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Fig. 5. — The Ha and H/3 activity fractions as a function of absolute vertical distance from the Galactic plane for M2-M7 dwarfs. For 
this comparison, only the stars with data good enough to measure H/3 were included. Error bars indicate the binomial errors in the Hq 
fractions — the H/3 uncertainties are similar. The H/3 activity fractions are in excellent agreement with the Ha data, indicating that H/3 
activity has the same duration in M dwarfs. 



Each panel of Figure 7 shows the median luminosity of 
a specific activity tracer (as compared to the bolometric 
luminosity of the stars) for active stars sensitive to that 
particular emission line. All of the activity tracers ex- 
hibit similar behavior, with the relative luminosity stay- 
ing constant for the early-type M dwarfs and then falling 
to lower values at or around a spectral type of M5. The 
median values and the upper and lower quartiles that are 
plotted in Figure 7 are also given in Table 4. The bluer 
emission lines produce smaller fractional luminosities as 
expected, and are consistent with previous Balmer decre- 
ment studies (Hawley & Pettersen 1991; Bochanski et al. 
2007b; West & Hawley 2008). W08 found evidence for 
a slight decrease in Lna/Lboi as a function of Galactic 
height, but the effect is small and not included in our 
analysis. 

3.3. Metallicity Sensitive Features 

Metallicity continues to be one of the most elusive 
physical quantities to measure in M dwarfs. Recent work 
using infrared spectroscopy and improved photometry for 
M dwarfs in wide binaries may soon lead to accurate 
metallicity determinations for large samples of M dwarfs 
(Woolf et al. 2009; Johnson & Apps 2009; Rojas-Ayala 
et al. 2010). However, the Lepine et al. (2007) C Pa- 
rameter is currently the best indicator of metallicity in 
optical M dwarf spectra. Using the DR7 M dwarf sam- 
ple, we examined how ( changes as a function of Galactic 
height. Figure 8 shows median ( values for Ml and M2 
dwarfs as a function of absolute height above the Galac- 
tic plane. These specific spectral types were chosen be- 



cause they span a large range of distance and contain a 
significant number of stars with ( uncertainty < 0.1. Fig- 
ure 8 shows a clear decrease in ( for both spectral types 
with increasing Galactic height. The discrepancies in the 
nearby bins reflect the scatter in the high-metallicity end 
of the C vs. metallicity relation (Woolf et al. 2009). The 
results shown in Figure 8 suggest a significant change in 
the median metallicity of M dwarfs (and by extension the 
Galaxy) as a function of age, with the oldest stars (that 
are farther from the Plane) having the lowest metallic- 
ity. Similar trends have been seen in both the radial and 
vertical directions from previous spectroscopic samples 
of higher mass stars (Nordstrom et al. 2004; Ivezic et al. 
2008). However, the higher density and longer main se- 
quence lifetimes of M dwarfs may provide a better probe 
of the metal content (and evolution) of the local Galactic 
disk. The ( values of ~0.6 in the most distant bins of 
Figure 8 have metallicities of [Fe/H]~-0.7 according to 
the wide binary analysis of Woolf et al. (2009).^^ We 
note that the photometric distances derived in our sam- 
ple are likely affected by metallicity and that the farthcist 
bins in Figure 8 (that have lower metallicities) are closer 
than shown due to fainter Mr magnitudes (Paper II). 
Regardless, this strong trend in metallicity demonstrates 
the utility of large spectroscopic samples of M dwarfs 
(such as the one presented in this paper) for understand- 
ing the chemical evolution of the Milky Way. However, 

19 The Woolf et al. (2009) [Fe/H] - C calibration was lim- 
ited to ~M0-M3 dwarfs and has considerable scatter at the high- 
metallicity end of the relation. 
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Fig. 6. — The Ha, H/3, H7, H<5 and Call K activity fra<;tions for M4 dwarfs as a function of absolute vertical distance from the Galactic 
plane. Each panel includes only the stars where the data are good enough to be sensitive to that emission line (as well as all of the 
redder emission lines). Uncertainties axe included for Ha but are similar for all of the emission lines. The activity fractions have the same 
magnitude and decreasing trend with Gala<;tic height for all lines. The a<;tivity fra<;tions for stars sensitive to H5 and Call K activity are 
slightly higher than for stars sensitive to H/3 and H7. 



TABLE 4 

Luminosity in Activity Trac;er Emission Lines 



Spectral 


^Ha/Lbol 


-^'H/s/Lbol 


-t'H7/Lbol 




iCall K/Lbol 


Type 


(xlO-4) 


(xlO-") 


(xlO-4) 


(xlO--*) 


(xlO-") 



MO 



L.27i;-? 0.85t;;;?« 0.50i»;i| 0.76i»:i« 0.33i°;?I 

Ml 1.51«;S 0.811";??, 0.50t;;;?« 0.39t»:i« 0.45t°;?i 

M2 1.42t°;« 0.77li;;ri OAit^]^ 0.32t»:«| 0.34t°;J| 

M3 1.36tO-g 0.67li',??, 0.35^;!-^? 0.26t'^;l° O.SOtVit 

M4 i.44i;;;g o.6iii;:i o.35i;;;i^ o.29t^^ o.28to;ii 

M^ nSl+"-2-"'' n 33+°-" (107+0.10 26+"-°^ Q 07+0.09 

IViO "-^^-0.22 '^••^■'-0.08 '^■■''-0.11 '-'■^"-0.07 '^■■''-0.15 

M6 0.571°;?° 0.221°;"^ 0.131";°? 0.08l»;°? 0.171°;°^ 

M7 0.l9l°;°i 0.06l°;»? 0.071°;°° 0.06l»;°l 0.05l°;°i 

M8 0.271°;;? 0.14l°;°« 

M9 o.2oi;;;°-^ 0.121°;;^, 



Note. — Bolometrically scaled luminosities were derived using stars 
that show activity in each emission line. Reported uncertainties are the 
25th and 75th quartile values. 



additional calibration is required to convert quantities 
such as C into precise metallicitics. 

Previous SDSS studies have shown that the g — r color 
of M dwarfs appears to correlate with metallicity (West 
et al. 2004; Lepine & Scholz 2008). Future large photo- 
metric surveys such as LSST (Ivezic et al. 2008) will rely 
heavily on the ability to classify objects based on their 
broadband colors alone. A large spectroscopic sample of 



M dwarfs, like the one presented here, can help corre- 
late photometric properties of M dwarfs to spectroscop- 
ically derived physical parameters. Figure 9 shows the 
g — r colors of M dwarfs as a function of ( for M1-M4 
dwarfs. Although there is a trend of decreasing ( with 
rc^ddcr g — r color, the change in g — r is different across 
spectral types. Therefore, g — r color alone is not a good 
tracer of M dwarf metallicity across the range of M dwarf 
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Fig. 7. — The ratio of line luminosity to bolometric luminosity as a function of spectral type for Ha, H/3, H7, H<5 and Ca II K. The higher 
order Balmer and Call K emission lines have progressively lower luminosity than the Ha lines. All of the emission lines show a drop in the 
luminosity ratio for late- type M dwarfs. 



spectral types. However, when we examined the ( val- 
ues as a function of both g — r and r — z (a proxy for 
spectral type; see Tabic 2), we found a strong trend of 
decreasing ( diagonal to the g — r, r — z axis and perpen- 
dicular to the stellar locus (Covey et al. 2007). Figure 
10 highlights a region of color space that can be used 
to probe the nictallicity content of M dwarfs and obtain 
low-mctallicity members of the Galactic halo from solely 
broadband photometry. Using a Lcvenbcrg-Marquardt 
minimization technique, we derived a two-dimensional 



fit that relates ( to the g — r and r — z colors: 



C = 1.04 - 0.98(5 - rf - 0.07(r - z)^ + 1.07{g 



-0.53(r -z) + 0.63(5 - r){r - z). 



(3) 



Equation 3 is valid over the ranges shown in Figure 10 
and has typical uncertainties of 10-20%. 

3.4. Nearby Stars 
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chemical evolution of the nearby Milky Way. 
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added catalog includes proper motions, RVs, photomet- 
ric matches to 2MASS, spectral classification, distances, 
calculated space velocities, activity-induced emission line 
measurements and molecular bandhead strengths. The 
DR7 catalog is available for download at the Vizier site 
or by contacting the corresponding author. Our anal- 
ysis of the visual spectral classification as compared to 
the automatic Hammer results reveals a slight ^ 0.4 sub- 
type systematic offset in the automatic Hammer spectral 
types for late-type M dwarfs but confirms that the auto- 
matic types are good to within the stated precision of 1 
spectral subtype. We present updated median colors for 
all M dwarf spectral types. 

We have also analyzed some of the bulk properties of 
the DR7 sample. The main results of our analysis are as 
follows: 

1. The magnetic activity fractions of low-mass stars 
as traced by H/3, H7, H5, and Call K, decrease as 
a function of Galactic height and agree with those 
previously traced using Ha. This confirms that 
the presence and duration of magnetic activity as 
traced by the higher order Balmer emission lines 
and Call K is similar to that of Ha, albeit at lower 
line luminosities and only measured in stars with 
sufficient SNR. 

2. The metallicity sensitive parameter ( decreases as 

a ftmction of Galactic height, confirming a decline 
in the metal content of distant M dwarfs (and the 
Galactic disk) as a function of age. 

3. As previously shown, the g — r color of M dwarfs 
correlates with metallicity (here parameterized by 
C). However, the g — r versus ( relation is spectral 
type dependent. Using both the g — r and r — z 
colors, we demonstrate that low-metallicity sub- 
dwarfs can be identified using photometry alone. 
This relation will be useful for source identifica- 
tion in upcoming large photometric surveys such 
as Pan-STARRS and LSST. 



Prompted by the discovery of new low-mass stars in 
the solar vicinity from other SDSS studies (e.g. Schmidt 
et al. 2010a), we inspected our catalog for M dwarfs that 
are possibly within 25 pc of the Sun. We found 21 nearby 
M dwarf candidates (8 of which were previously uniden- 
tified) based on the distances derived using the Adr, r — z 
photometric parallax relation of Bochanski et al. (2010). 
Three of the stars were in the W08 SDSS DR5 sample 
but were not identified there as potential solar neighbors. 
Table 5 gives the positions, spectral types and distance 
estimates for all 21 of the nearby candidates. The clos- 
est candidate is SDSS1410+1846, which has an estimated 
distance of 14.8 pc. Since roughly one out of every 1000 
M dwarfs in SDSS has a spectrum, there is a strong pos- 
sibility that many more nearby stars await discovery and 
will be cataloged in future studies. 

4. CONCLUSIONS 

We have presented the SDSS DR7 M dwarf spectro- 
scopic catalog, which consists of more than 70,000 visu- 
ally confirmed M dwarfs and represents the largest spec- 
troscopic sample of M dwarfs ever assembled. Our value- 



4. The DR7 M dwarf sample contains several previ- 
ously unidentified M dwarfs that are likely within 
25 pc of the Sun, including one that is possibly 
closer than 15 pc. The sparse spectroscopic cover- 
age of low-mass stars in SDSS suggests that there 
arc numerous M dwarf solar neighbors that will be 
identified in future studies. The advent of large, 
multi-epoch, deep surveys will be particularly use- 
ful for completing the nearby M dwarf census. 

Many additional studies will make use of our M dwarf 
catalog. Those already underway include a statistical 
parallax analysis of absolute magnitude variations in the 
M dwarf population (Paper II), a detailed examination 
of M dwarf kinematics and the motions of the Galactic 
thin and thick disks (Pineda et al. 2011, Paper HI), and 
an investigation of the content and distribution of dust 
in the local Galaxy (Jones et al. 2011). In addition to 
the current and unforeseen science that will be accom- 
plished with our sample, we anticipate that our M dwarf 
catalog will be used to select and classify M dwarfs in 
several upcoming large surveys. Over 1 billion M dwarfs 
will be observed and cataloged in the new wave of large 
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Fig. 10. — The g — r versus r — z colors of M dwarfs in the DR7 spectroscopic sample. The bins have been color coded according to 
their f values. Metallicity decreases along roughly diagonal lines in the g — r, r — z color space and perpendicular to the stellar locus 
(Covey et al. 2007). This relation should be useful for classifying M dwarfs and identifying low-metallicity subdwarfs in upcoming large 
photometric surveys. 

photometric surveys coming online in the next decade. 
We hope that our large M dwarf sample presented in 
this paper will provide a useful tool for correlating the 
spectroscopic attributes of low-mass stars with their pho- 
tometric properties both in single exposures and in the 
time domain. 
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